Density of states of disordered systems with a finite correlation length 
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We consider a semiclassical formulation for the density of states (DOS) of disordered systems in 
any dimension. We show that this formulation becomes very accurate when the correlation length 
of the disorder potential is large. The disorder potential does not need to be smooth and is not 
limited to the perturbative regime, where the disorder is small. The DOS is expressed in terms 
of a convolution of the disorder distribution function and the non-disordered DOS. We apply this 
formalism to evaluate the broadening of Landau levels and to calculate the specific heat in disordered 
systems. 
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Strong disorder plays an important role, particularly 
in amorphous systems such as semiconductors or glasses. 
In electronic systems and in waves in general, localiza- 
tion due to disorder is common and has been studied 
for many years Q, |^ HI . It was shown that localization 
can strongly affect transport properties, such as inducing 
metal-insulator transitions. By contrast, thermodynam- 
ical properties like the specific heat only depends on the 
DOS and not on the physics of localization. 

Evaluating the DOS for disordered systems is usually 
very difficult and is often done with weak disorder ex- 
pansion techniques. However, some exact results exist 
and include the smoothing out of band edges or singular- 
ities 1^ and the existence of Litshitz tails In general, 
the DOS decays exponentially close to the band edge 
for most disorder distributions 0]. However, there are 
no general results valid for the entire DOS in the case 
of strong disorder and particularly in dimensions larger 
than one. Furthermore, most of the results are based on 
the assumption of an uncorrelated disorder potential. In- 
deed, when a non-zero correlation length ^ is introduced 
in the disorder potential the situation is not well under- 
stood. This is nonetheless important for the understand- 
ing of many physical systems, since the correlation length 
is often non-zero. For example, in quenched disordered 
systems like amorphous metals the nature of the disorder 
is largely dependent on the cooling rate and subsequent 
annealing procedure 0- In amorphous semiconductors 
the correlation length depends on deposition conditions 
[3, whereas in GaAs/AlGaAs heterostructures the cor- 
relation length can be modified by bias cooling 0. 

In this letter we reformulate the semiclassical deriva- 
tion of the DOS and show that it is a very good approx- 
imation to the exact result in the limit where the cor- 
relation length ^ of the disorder potential is large. Our 
approach is based on a semiclassical formulation, which 
originated with Kane |lfij | and was further developed by 
Ziman • More recently, a semiclassical approach was 
used in the context of amorphous semiconductors [l^ 



and in ID Here we introduce a more general situa- 
tion, where we consider systems, where the volume t{E) 
of energy E in the phase-space can be partitioned for any 
given realization like: 



t{E)=Y^t({E-V,), 



(1) 



where r/ {E — Vi) is a volume of energy E — Vi oi the 
partition i of the total volume t{E). Here r/ corresponds 
to the volume free of disorder, since in this partition Vi is 
assumed approximately constant for a given realization 
of disorder. 

While the assumption in |^ is quite general and 
can be applied to a number of systems, we will focus 
our attention on the case of a particle with energy E 
in a random potential. The potential is constant by 
pieces Vi over a constant size d, which represents the 
correlation length, ^, of the disorder potential. The 
Schrodinger equation in any dimension is then simply 
Eip = Hfip -f ip Vi{x — Qi), where Hf is the Hamilto- 
nian of the disorder free system and the functions Vi are 
random independent functions with a given probability 
distribution centered at dl. The potential energy Vi is 
approximately constant for a given realization. In this 
case, the total volume t{E) can be written as a function 
of the sub- volume like in eq.Q). 

Semiclassically, the total number of states N{E) is then 
proportional to the volume t(E) in phase space. Namely, 
N{E) = T(E)/h This is our main starting point to 
obtain an expression for the DOS of disordered systems. 
Further below, we will then analyze under which con- 
ditions (0 holds. We assume that {Vi} is an ensemble 
of random independent variables with known probability 
distribution Ptot ~ Y[Pi^i)- The disorder average (•) of 
(HJ is given by (t{E)) = J^ii'^I ~ ^0) and since the 
distribution of the random independent variables Vi is 
the same, we have {t{E)) = / dVPiV) |r/ {E~V)^ 
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and {t{E)) = JdVP{V)\J2^'r^iE~V)y Because 

the summation t[ {E — V) corresponds to the vol- 
ume of the disorder-free system, namely, t^{E — V) — 
~ the expression for the integrated DOS 

N{E) = T{E)/h becomes 

{N{E)) = J dVP{V)Nf [E - V) (2) 

and for the DOS g{E) = dN/dE we have 

{g{E))^ J dVP{V)9^{E-V). (3) 

The disorder-free case corresponds to Dirac's distribu- 
tion P{V) = 5{V), which leads to {g(E)) = gf (E) as 
expected. This expression can also be obtained by us- 
ing the summation of local density of states following 
Kane's original work [lOj . In what follows we analyze 
under what conditions expressions (|2I3|I can be used and 
when they are accurate. We will start by considering 
the Anderson model Q, i.e., a discretized version of the 
Schrodinger equation, which is X^j^yV'j = ^ Vi)'ipi, 
where Tij = — 1 for nearest neighbors and otherwise 
and Vi is our random potential. The DOS of the non- 
disordered cubic case is then simply given by 

D 

g{E)^ / 5(£; + 2 Vcos(fc,))rf^k, (4) 

where D is the dimension. While evaluating the integral 
Q gives no simple expression for g{E) for higher dimen- 
sions, g{E) is nonetheless a symmetric function of E and 
the DOS near the band edges Eq = ±2D scales simply 



as g{E) ^ \E — i?o|^^~^^^^ like in the continuous case. 
When introducing disorder, the DOS broadens. 

We first considered the case of uncorrelated disorder, 
i.e., (ViVj) ~ Si J, and where the l^'s are uniformly dis- 
tributed between —a and ct. For this case (D=3) there ex- 
ists no analytical result for the entire shape of the DOS. 
Hence, we evaluated the DOS numerically and show the 
result in figure ^ for different strengths of disorder. In 
the same figure we also show the DOS evaluated with 
Clearly the agreement is not good. This is not too 
surprising and will be discussed below. 

The situation changes quite dramatically, when we in- 
clude the condition for the derivation of our semiclas- 
sical expression Indeed, when we impose that the 
potential is constant by pieces, we obtain a much better 
agreement between expression Q and the exact diago- 
nalization (see figureEJ. In figureElwe show the 3D DOS 
for the case where Vi is constant over a fixed interval, 
ranging from 1 to 4. The larger the interval the better 
the agreement. It is important to note that increasing 
the interval d does not reduce the strength of the disor- 
der (V^), but it only increases the disorder correlation 
length, where £, = d. While the overall quality of the fit 
is quite remarkable as soon as d exceeds 2, the band tails, 
however, are only accurate for \E — Emin\ > SE, where 
Emin corresponds to the band edge of expression SE 
can be estimated from the cut-off of the wavelength due 
to d, hence 6E ~ E{kmin + 2TT/d) — E{kmin), where E{k) 
is the dispersion relation of the disorder free system. We 
checked this by evaluating the ID case, which is shown in 
the inset of figure |51 where the accuracy of the band tail 
DOS is improving drastically for d increasing from 4 to 
50. It is also interesting to note that the sharp structures 
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FIG. 1: Dotted lines: average density of states calculated for 
a 3D system of size 16^ and averaged over 20 configurations 
for the Anderson model with a uniform distribution of width 
a = 4,8, 16. The bandwidth of the non-disordered case is 12 
in these units. The solid lines are the expressions obtained 
from (Ol. 



FIG. 2: Average density of states calculated for a system of 
size 16^ and averaged over 30 configurations for the Anderson 
tight binding model Q with a uniform distribution of width 
CT = 16 and d = 1,2, A. The black line is the expression 
obtained from l^l. In the inset the ID density of states is 
shown for d=4, 16, 50 and the black line is obtained from iPJ. 



3 



close to the band center are accurately described by Q 
for d large enough. 

The disorder potential considered above, while long 
range in nature, is different from the often consid- 
ered "smooth" disorder case. For comparison, we 
considered a smooth potential of the form: Wi = 
(dv^)"'/'Ej^je-('-J')'/'*', where (W^) = {V^) and 
the correlation length (FWHM) is given by ^ = 2.35d. 
The distribution function for this potential at large d 
is simply given by P(W) = v^27r(V^2^e-^'/(2<^'». We 
evaluated numerically the exact DOS in ID and com- 
pared it to expression Q for different values of d using 
the exact numerical distribution function. The results 
are shown in figure 13 Here again, expression (jS)) is in ex- 
cellent agreement with the exact result when d is large. 

For smooth disorder this is not too surprising, since the 
semiclassical approximation is accurate for a very smooth 
potential. Indeed, a recent perturbative calculation in 
ID shows that the semiclassical approximation is exact 
in the limit, where {(VV)^) / (V^)^/^ < 1 [Ti. They 
also obtained the first order correction of the DOS to the 
semiclassical result. In this letter, however, we push the 
argument much further and show that for a disordered 
system the semiclassical expression Q is accurate in the 
limit where the correlation length is large, independently 
of the smoothness and dimension of the system under 
consideration. Therefore, we have an expression, which 
provides a simple form for the DOS of disordered sys- 
tems with a large correlation length. This is the main 
result of this letter. In localization physics, disorder cor- 
relations can also have a dramatic effect on the extend 
of the wave functions and even lead to the existence of 
extended states in ID |0| and 2D 0|. Without correla- 
tions, Lloyd showed that for P{V) a Lorentz distribution, 
the exact DOS is also given by expression Q 0|. 




FIG. 3: Exact average DOS in ID for a smooth disorder po- 
tential with characteristic length d — 1 and d — 16 and the 
DOS obtained from . The distribution of the random po- 
tential for d = 16 is shown in the inset. 



We now turn to some important consequences of ex- 
pression Q. We start with the singularities in the 
disorder-free system, which will be smoothed out if the 
disorder distribution is continuous. An interesting ex- 
ample of this is the situation of a perpendicular mag- 
netic field in two dimensions, where the spectrum is dis- 
crete and the density of states is given by a sum of S 
functions centered at the Landau levels, i.e., g{E) ~ 
J2n^i-^ — (n + l/2)hajc), where iVc is the cyclotron fre- 
quency. Using expression Q to obtain the average DOS 
in the presence of disorder immediately gives a broad- 
ening of the Landau levels only related to the disorder 
strength and the disorder distribution. Hence, 

V2)?ic.c) (5) 

n 

It is interesting to note that this Landau level broadening 
does not depend on the Landau level index. While the 
exact shape of the Landau level broadening is still under 
debate, experiments on GaAs/AlGaAs heterostructures, 
which have long range disorder, indicate that the broad- 
ening is magnetic field dependent at low fields, where the 
Landau level overlap is significant. At higher fields the 
broadening becomes independent of the field and of the 
Landau level index and the shape is found to be best fit- 
ted by a Lorentzian Hence, expression Q is consis- 
tent with experimental findings for large fields in systems 
with long range disorder. 

Another interesting consequence of expression Q is 
the widening of the bands. Indeed, if P{V) has a finite 
support, namely, —a < V < +a, then the bands will be 
extended by a on each side. This can lead to the decrease 
of gaps in disordered materials. The band edges will ac- 
quire tails. Expression ^ leads to tails, which depend 
on the disorder distribution. For continuous systems the 
disorder free band edge is given by g{E) - in D = 3. 
Introducing disorder will modify this dependence. It is 
quite straightforward to see that, for a square distribu- 
tion P{V) = l/2cr for -a < V < +(7 and P{V) = 
otherwise, we obtain in any dimension 

(g{E)) = ^ {Nf{E + Nf{E ~ a)] , (6) 

where is the integrated DOS of the disorder-free sys- 
tem. Therefore, {g{E)) ^ {E + tr)^/^ at the band edge 
and the band tail is not exponential for this distribu- 
tion. Numerically, we see that this is indeed the case for 
the correlated case (c? = 50) shown in the inset of figure 
12 For other disorder distributions, P{V), such as Gaus- 
sian or Poisson, we do recover an exponential band tail. 
Fundamentally, the disorder distribution will govern the 
dependence of the band tails. In the high disorder limit 
E/a < 1, {g{E)) ~ (ct) is independent on energy. 

We now turn to the specific heat of disordered systems. 
The internal energy Up and the particle number Np 
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for fermions is given by Up = J ^(e-J^ikt^^^ {9{E))dE 
and Np — J ^(e-„))kt^i {9{E)) dE. As above we con- 
sider the disorder-free DOS g{E) = a-ji?~2— ^ where 
a is a constant. In the high temperature Hmit (Boltz- 
mann distribution) and using © leads to a specific heat 

CboU. = ^ or 



do recover that the specific heat is smaller in the dis- 
ordered case than in the ordered case. In the opposite 
limit of small relative disorder the situation is reversed 
because the disorder-free DOS in ID decreases with en- 
ergy, leading to a comparatively larger specific heat for 
the disordered case. 



CboHz — K ■ 



D 



KTsmha/KT 



(7) 



For the case Z? = 1 we recover the results found in ref- 
erence 01 . For all dimensions, the limit cr ^ gives 
CboHz = KD/2 in agreement with the Dunlap-Petit law. 
The case a ^ oo leads to Csoitz = K{D + 2)/2, which 
represents an enhancement of the specific heat due to 
disorder. 

In the opposite limit of low temperatures, where we can 
use the Sommerfeld expansion, the difference between the 
specific heat of the disordered system and the ordered 
system (for the same electronic density) is given by 



In conclusion, we have shown that the average DOS of 
disordered systems can be derived in a semiclassical way 
to represent an accurate DOS for a disorder potential 
when the correlation length is large. We further, evalu- 
ated a few physical quantities such as the broadening of 
the Landau levels, the characterization of band edges and 
tails and finally, discussed the electronic specific heat in 
these systems. 



Cn — Cv 



3 12 



(2 - D)D, 



(8) 



where E^ is the Fermi level of the disordered system and 
we further assumed that a <^ Ejj. Clearly this expres- 
sion shows that the disorder decreases the specific heat 
in 3D but enhances it in ID. The 3D case is simple to 
understand, since at low temperatures the specific heat 
is mainly determined by the tail of the DOS. Hence, for a 
given energy bandwidth, the DOS of the disordered case 
is smaller than the DOS of the ordered case, which leads 
to a smaller specific heat in the disordered case. The ID 
result is more surprising, but can be understood by look- 
ing at all relative disorder strengths presented in figure 
01 Indeed, at high disorder strength, (i.e. Ed small), 
where the low energy DOS is dominated by the tails, we 
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FIG. 4: Difference between the specific heat of the disordered 
system and the ordered system in units of ^^^l "^ Q as a func- 
tion of the Fermi level of the disordered system. 



M.H. would like to acknowledge support from RQMP, 
FQRNT and NSERC. 



[9 

[lo; 
[11 

[12; 

[is: 

[14 

[is: 



[16: 

[17: 
[is: 



P.W. Anderson, Phys. Rev. 109, 1492 (1958). 

P.A. Lee and T.V. Ramakrishnan, Rev. Mod. Phys. 57, 

287 (1985). 

CM. Soukoulis and E.N. Economou, Wave Random Me- 
dia 9, 255 (1999). 

F.J. Wegner, Phys. Rev. B 19, 783 (1979). 

I.M. Lifshitz, Sov. Phys. Usp. 7, 549 (1965). 

L. Pastur and A Figotin, Spectra of Random and Almost 

periodic operators. Springer Verlag, (1992). 

S.R. Elliott, Physics of amorphous materials, John Wiley 

& Sons (1989). 

D. Quicker and J. Kakalios, Phys. Rev. B 60, 2449 
(1999). 

PT. Coleridge, Semicond. Sci. Technol. 12, 22 (1997). 

E. O. Kane, Phys. Rev. 131, 79 (1963). 

J. M. Ziman, Models of Disorder (Cambridge University 
Press, 1979). 

S.K. O'Leary, S. Zukotynski and J.M. Perz, Phys. Rev. 
B 51, 4143 (1995). 

J. C. Flores, Phys. Stat. Sol. (b) 239, 168 (2003). 

D.N. Quang and N.H. Tung, Phys. Rev. B 60, (1999) 

13648. 

J.C. Flores, J. Phys. Condens. Matter 1, 8471 (1989) and 
D.H. Dunlap, H.-L. Wu, and P.W. Phillips, Phys. Rev. 
Lett. 65, 88 (1990). 

M. Hilke, J. Phys. A: Math. Gen. 27, 4773 (1994) and 

M. Hilke, Phys. Rev. Lett. 91, 226403 (2003). 

P. Lloyd, J. of Phys. C 2, 1717 (1969). 

A. Potts et al, J. of Phys.: Condens. Matter 8, 

5189(1996). 



